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Abstract 

The various equations at the surfaces and triple contact lines of a deformable body are 
obtained from a variational condition, by applying Green's formula in the whole space and 
on the Riemannian surfaces. The surface equations are similar to the Cauchy's equations 
for the volume, but involve a special definition of the 'divergence' (tensorial product of the 
covariant derivatives on the surface and the whole space). The normal component of the 
divergence equation generalizes the Laplace's equation for a fluid-fluid interface. Assuming 
that Green's formula remains valid at the contact line (despite the singularity), two equa- 
tions are obtained at this line. The first one expresses that the fluid-fluid surface tension is 
equilibrated by the two surface stresses (and not by the volume stresses of the body) and 
suggests a finite displacement at this line (contrary to the infinite-displacement solution of 
classical elasticity, in which the surface properties are not taken into account). The sec- 
ond equation represents a strong modification of Young's capillary equation. The validity 
of Green's formula and the existence of a finite-displacement solution are justified with an 
explicit example of finite-displacement solution in the simple case of a half-space clastic solid 
bounded by a plane. The solution satisfies the contact line equations and its elastic energy is 
finite (whereas it is infinite for the classical elastic solution). The strain tensor components 
generally have different limits when approaching the contact line under different directions. 
Although Green's formula cannot be directly applied, because the stress tensor components 
do not belong to the Sobolev space i? 1 (V), it is shown that this formula remains valid. As a 
consequence, there is no contribution of the volume stresses at the contact line. The validity 
of Green's formula plays a central role in the theory. 
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1 Introduction 



Surface properties of deformable bodies have been continually studied since the early work of 
Gibbs [6] until recent mechanical or thermodynamic studies, e.g. [H[18j|21][Ml[25]. They have 
many applications, e.g. in adhesion, coating and nanosciences (since small and thin objects are 
deformable and have a high surface/volume ratio). A previous paper [18] was devoted to the 
physical basis of the theory: application of the equilibrium criterion of Gibbs; introduction of the 
new concept of 'ideal transformation', i.e., the homogeneous extrapolation of the deformation, 
in the interface film, up to the dividing surface; determination of the thermodynamic variables 
of state of a surface; definition of the surface stress tensor; surface and line equations. Moreover, 
for an elastic solid, it is known that classical elasticity predicts a singularity with an infinite 
displacement (and an infinite elastic energy) at a solid-fluid-nuid triple contact line, owing to 
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the fluid-fluid surface tension which is a force concentrated on this line |22tl23|. Although 
some authors tried to overcome this problem (by introducing some fluid-fluid interface thickness 
|12U2U| . some cut-off radius near the contact line [22, 23j or some new elastic force at this line 1 1 '■)[ ) . 
this situation makes very difficult to write any equilibrium equation at the contact line. 

The present paper concerns the mathematical foundation of the theory. A sketch of the proof 
of the surface and contact line equations is presented (no proof was given in the previous physical 
paper [IB]), which shows (i) the importance of the validity of Green's formula at the contact 
line (despite the singularity) and (ii) owing to the surface properties, the probable existence 
of a finite-displacement solution (consequence of the line equations, based themselves on the 
assumption of the validity of Green's formula). These two points are justified with an explicit 
example of finite-displacement solution, in the simple case of a half-space solid, bounded by 
a plane, and subjected to a normal force concentrated on a straight line of its surface. This 
solution also shows that the elastic energy is finite and that Green's formula remains valid at 
the contact line. 

2 Surface and Contact Line Equations 

For a general deformable body b in contact with various immiscible fluids f, f',... (with no mass 
exchange between the body and the fluids), the mechanical equilibrium condition relative to the 
body, including its body-fluid surfaces (bf, bf',...) and its body-fluid-fluid triple contact lines 
(bff ',...), may be written as 



(: means double contraction; see [THJ for the physical basis of the theory), in which 5 is an 
arbitrary variation such that, on the closed surface X which bounds the system, the points of 
the body and the points of the body-fluid-fluid lines remain fixed. In this expression, tt is the 
Piola-Kirchhoff stress tensor (i.e., the Lagrangian form, relative to a reference state of the body) 
at equilibrium, e the Green-Lagrange strain tensor (also relative to this reference state), dv and 
dvo are respectively the volume measures in the present state and in the reference state, p is 
the mass per unit volume, g the (constant) gravity vector field, 5x the displacement of a point 
of the body, p the fluid pressure, n the unit vector normal to the bf surface, oriented from f 
to b, da and dao are respectively the area measures in the present state and in the reference 
state, p s is the mass per unit area (excess on the dividing surface Sbf defined by the condition: 
no excess of mass of the constituent of the body), ir s the (Lagrangian) surface stress tensor at 
equilibrium, defined in [THJ, e s the (Lagrangian) surface strain tensor, defined in Appendix |A"1 
7fP the fluid-fluid surface tension, v^i the unit vector normal to the bff' line, tangent to the ff' 
surface, and oriented from the line to the interior of ff', 5X the (vector) displacement of the 
bff' line, perpendicular to the line (in the present state), dl and dlo are respectively the length 
measures in the present state and in the reference state, 70 is the surface grand potential (excess 
on the dividing surface), per unit area in the reference state, and 5Xq the (scalar) displacement 
of the bff' line, measured in the reference state, perpendicular to that line in the reference state, 
and positively considered from bf to bf (see also Fig. [21 below). 

This variational equilibrium condition leads to various equations at the surfaces and the 
triple contact lines of the body. Since the preceding paper [THJ was devoted to the physical 
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aspects of the theory, these equations were only written without proof. In this section, a sketch 
of this proof is presented, which shows the importance of the validity of Green's formula to 
obtain the contact line equations (despite the line singularity). These equations then suggest 
the existence of a finite-displacement solution. 

In order to only have quantities or variables (such as points, forces, etc.) which refer to the 
present equilibrium state in these equations, we first need to transform all the Lagrangian terms 
in (|2.ip (i.e., those related to the 'undeformed' reference state) into Eulerian forms (i.e., related 
to the deformed present state). It is well known that the Eulerian forms of the above (volume) 
stress and strain tensors, ir and 6e, are the Cauchy stress tensor a and the infinitesimal strain 
tensor 5e defined below in the next paragraph (see e.g. [14], tome I, annexe II). Note that e 
measures the strain between the 'undeformed' reference state and the deformed present state, 
so that its components may have arbitrary values, since large strains may occur in highly 
deformable bodies (even when subjected to capillary forces or surface stresses). Note also that 
5e, which measures the infinitesimal strain between the present state and its varied state (i.e., 
after the variation 8), is not the variation of some strain tensor, but it is related to the variation 
5e of the Lagrangian tensor e by 

5e = $Q-5e- $ , 

where 3>o is the deformation gradient between the reference state and the present state and <I>g 
its adjoint. Classically, the work of deformation of a volume element (first term of (|2.1|) ) may 
be written in the Eulerian form 

7r : 5edvo = a : 5e dv (2-2) 

(see [14| . ibid.), i.e., with arbitrary Cartesian coordinates in the three-dimensional space E 

ir 13 Seij dvo = a %3 5eij dv, 

where Latin indices i, j, k,... belong to {1,2,3} and summation is performed over repeated 
indices. In a similar way, these concepts are extended to the surfaces in Appendix [Aj where 
the Lagrangian surface strain tensor e s , the Eulerian infinitesimal surface strain tensor Se s and 
the Eulerian surface stress tensor a s are defined. The work of deformation of a surface element 
(fifth term of (|2.ip ) may then be written in the Eulerian form (|A.12p . 

Let us first consider the simple case of a bounding surface £ which only encloses one fluid f 
and the body b. The equilibrium condition (|2.ip may then be written as 



/ a : be dv — pg ■ w dv 
Jv Jv 

— / pn ■ w da — / p s g ■ w da + / er s : Se s da = 0, 
Js Js Js 

where w = Sx, V is the bounded open set of E occupied by the part of the body enclosed in S, 
and S the bounded part of Sbf enclosed in S. Since Se = s ((Dm)* + Dw) and 

f tr(a-5e)dv = [ tr( — ^-^--Dw)dv 
Jv Jv 2 

= / tr(cr* • Dw) dv + / tr( - ° -Dw)dv, 
Jv Jv 2 

by application of Green's formula (with w = on S) 

/ tr(cr* • Dw) dv = — div(cr*) ■ wdv — / (a* ■ w) ■ nda 
Jv Jv Js 

div(<r*) ■ wdv — l (a ■ n) ■ wda (2-3) 
v Js 
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(if the components of a and w belong to C 1 (V); e.g. [2J, Sec. 3.2.1), this leads to the classical 
Cauchy's equations for the body 

diva + pg = (2.4) 
a* = a, (2.5) 

where diva is the vector associated to the linear form div(<7*), and the remaining condition for 
the surface 

— / (a ■ n) ■ w da — / pn ■ w da — / p s g ■ w da + / a s : 5e s da = 0, (2-6) 
Js Js Js Js 

for any variation such that w = on the closed curve T = Sbf H S which bounds S. Note that, if 
volume moments M dv were present, the new term - / v §tr(M* • Bw) cfa0 would appear in the 
equilibrium condition, and (|2.5p would become a — a* + M = 0. 

By applying Green's formula on the Riemannian manifold [3j S to the last term of the 
equilibrium condition, we then obtain 



(a ■ n) ■ wda — / pn ■ wda — I p s g ■ wda 
Js Js 

- [ div(5T) • w da + / ti C* ~ a * ■ i* ■ rp) da = (2.7) 
Js Js 2 

(see notations in Appendix where al = i ■ a s and the special divergence of al* is based on 
the tensorial product of the covariant derivative on T(Sbf) and the usual derivative on Sbf x E* 
(see Appendix [B]). This leads to the following equations for the surface 

div <f s + p s g + o ■ n + p n = (2-8) 
at = a s , (2.9) 

where div<f s is the vector associated to the linear form div(S\J*), i.e., according to t|B.4j) . 

dp{of dax') + ofT dax* + Ps g l + a* n 3 + pn' = 

°t a = of (2-10) 

(in the two last equations, a B and a as contravariant tensors, by raising the covariant index to 
the second place). Note the similarity of these equations with the classical Cauchy's ones (|2.4p - 
(|2.5p for the volume. Note also that (|2.9p might be different if surface moments were present 
(as for the volume stress: see the comment after (|2.6p ). The above Eq. (|2.8p has a tangential 
component 

div(j s + p s g t + (a ■ n) t = (2.11) 

(div a s being the usual surface divergence; the subscript t indicates the vector component tangent 
to Sbf) and a normal component 

l n : <7 S + p s g n + o nn +p = 0, (2.12) 

where l n = l-n (l n>a p = l l aj3 nf, I is the second vectorial fundamental form on Sbf), g n = 9' n an d 
<T nn = (<7-n) -n (see Appendix[B]). At any point x € Sbf, the eigenvalues of l n (as endomorphism) 
are the principal curvatures, and of Sbf ( [5], (20.14.2); a curvature being positive when 
its centre is on the side of n). Thus, if cr s is isotropic, i.e., a s = <f s I (eigenvalue d s and I the 



1 M is the endomorphism defined by (M ■ x) ■ y — [M, x, y], for any vectors x and y £ E; it satisfies M* = — M. 
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Fluid V 




Body b 



Figure 1: The bounding surface E encloses the parts V, S, S' and L of, respectively, the body b, the 
surface bf , the surface bf and the triple contact line bff' (the part of E in contact with the fluids and 
the surface ff' are not represented). 

identity), l n : a s = <f s tr(/ n ) = cf s (^" + t|j)> which shows that (|2.12)) generalizes the classical 
Laplace's equation for a fluid-fluid interface. 

Similar surface equations were obtained from a balance of momentum or equilibrium of 
forces [Tl[9lll51[2"4"]. a virtual power method [3], a thermodynamical approach [T]lllj. or an energy 
minimization [8ll25j. In these works, the existence of a surface stress tensor was often assumed, 
deduced from a given surface traction field [7J, or defined for elastic solids from a given set of 
thermodynamical or mechanical variables of state of the surface [TJ[8l[TTJ[25]- Note that our 
thermodynamical method [18], valid for any deformable body (not restricted to elasticity) and 
based on the general equilibrium criterion of Gibbs, leads to the determination of the 'local' 
thermodynamical variables of state of the surface, the definition of the surface stress tensor and 
the above equations. 

Let us now apply the equilibrium condition (|2.ip with a bounding surface £ which encloses 
two fluids, f and f, and the body b, in contact. V denotes the bounded open set of E occupied 
by the part of the body enclosed in E, S the bounded part of Sbf enclosed in E, S' the bounded 
part of Sbf enclosed in E, and L the part of the bff' triple contact line enclosed in E (Fig. [Q). We 
follow the same method as above, but Green's formula (|2.3p cannot be directly applied on V (S 
being here replaced with S U S'), owing to the singularity at the contact line. If b is a deformable 
solid subjected to a force concentrated on a line of its surface (here, the fluid-fluid surface 
tension 75/ applied on the contact line), then classical elasticity predicts a singularity with an 
infinite displacement at this line, together with an infinite value of the elastic energy |22^23j. 
Nevertheless, we shall see, in this paper, that the introduction of the surface properties leads to 
a solution with a finite displacement at the contact line and a finite value of the elastic energy. 
In the example of finite-displacement solution presented in the next section, the singularity at 
the contact line involves components of a which do not belong to H (V). Although Green's 
formula cannot be directly applied in this case (we would need that components of both a and 
w belong to H l (V); e.g. [2], Sec. 4.3.3), we show in Sec. 13.61 that this formula remains valid. We 
may thus assume that Green's formula is valid and, following the above method, the remaining 
condition for the surfaces and the line becomes (with the help of f|2.4[) - (|2. 5j) ) 



(a ■ n) ■ w da — pn ■ w da — I p s g ■ w da + a s : 5e s da 
sus' isus' jsus' Jsus' 



/ 7 ff / u s > -5Xdl+ / (7o,bf - 7o,bf) SX dl = 
Jl JLo 



'L 

(Lq is the position of L in the reference state), for any variation such that w = on the curves 
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Body b 




Reference state 



Present state and varied state 



Figure 2: Displacement 8Xq of the contact line bff' with respect to the body, in the reference state, and 
displacement SX of this line in the space, between the present state and its varied state, due to both 
the displacement of the line with respect to the body and the displacement of the material points of the 
body (see text). 

r = Sbf H £ and V = Sbp H £ which bound S U S', and the two points of L which belong to S 
remain fixed. Note that there is no singularity at Lo in the reference state of the body (which 
is, e.g., a state of the body before its contact with the fluids f and f ). The application of (jB.ip 
and Green's formula (jB.6[) to the two terms fa cr s : 5e s da + / g/ cr s : 5e s da leads to the two new 
terms — J L ((Ts,bf • ^bf) • ^bf dl — / L (o" s ,bf • ^bf) • ww dl (subscripts bf and bf respectively denote 
the sides of S = bf and S' = bf; thus, is the unit vector tangent to Sbf, normal to L and 
directed to the inside of Sbf; similarly, for with respect to Sbf; see Fig. [U w = on T and 
T', but not on the bf and bf sides of L) and then, with the help of ()2.8j) - f)2.9|) on S and S', to 
the remaining line condition 



for any variation such that the two points of L which belong to £ remain fixed (both in the 
space and with respect to the body b). Since the displacement SX of the contact line in the 
space is due to both the displacement of the corresponding material points of the body (w^f and 
Wbf on the bf and bf sides, respectively) and the displacement of the line with respect to the 
body (SXtf = 0o,bf ■ SXq and SX^i = 0o,bf • ^Xq on the bf and bf sides, respectively; here, 5X 
and 5Xq are considered as vectors, not necessarily normal to L and Lo, respectively; cfio defined 
in Appendix[A|), i.e., SX = u>bf + SX^ = ww + SX^i (see Fig. [2]), this condition becomes 



(the last term of the condition being written in Eulerian form, using 7 da = 70 dao for bf and 
bf), which leads to two equilibrium equations at the contact line (as in the case of the thin 
plate (IBI IT7]). The first one 



corresponds to a contact line fixed on the body (5Xq = 0, hence dX^f = SX^i = 0) and expresses 
the equilibrium of the two surface stresses and the fluid-fluid surface tension. Note that this 






L 



(7bf ^bf • SXu + 7bf vm> ■ 5X hi >) dl = 



C"s,bf ■ ^bf + cr s ,bf ■ Vbf> + 7ff' Vft> = 



(2.13) 
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equation suggests that a finite displacement occurs at the contact line (in the next section, an 
explicit example of finite-displacement solution will be presented). Some authors [13J proposed 
the presence at the contact line of a force originating from the volume stresses a in the body. 
The preceding equation shows that there is no such volume stress contribution. This is a 
consequence of the validity of Green's formula, as mentioned above, and will be illustrated in 
the next section (Subsec. I3.6p . With the help of (j2. 13fl . the above line condition gives the second 
equation (according to 5Xy,f = (fo,bf • $Xq an d SX^f = </>o,bf " <^o) 

^o,bf • (°"s,bf - 7bf I) ■ v\A + (f>o,w • (^s.bf - Tbf I') • fbf = (2.14) 

(/ and /' are the identity mappings on T a: (Sbf) and T x (Sbf)i respectively), which corresponds 
to a line moving on the body (SXq ^ 0), i.e., with <f> r = (f>o,bf' ' <^obf ( w hich is the 'relative 
deformation gradient' of the bf side with respect to the bf side; this concept was defined in |19j ; 
note that </> r does not depend on the reference state [18]). 

(er s ,bf - Tbf I) ■ fbf + (/)* ■ (<7 Si bF - Tbf /') ■ vw> = 0. (2.15) 

Applying t* to this equation (where r is a unit vector tangent to the contact line at x; thus, 
<j) x ■ r = t) gives the same equation as the tangential component (along r) of (|2.13|) . Applying 
^bf to (I2TT51) leads to 

&bf,vv — Tbf — {&W,vv — Tbf) &r,vv + <?W',tv «r,Tf = 0, (2-16) 

where (crbf,w> &m,tv) are the components of cr Sj bf ■ ^bf in the basis (vb{,r), similarly for cr Si bf 
with the basis (z4jp,t), and (a TtUU , a TjTU ) the components of 4> T ■ fbf in the basis (— vm',t) (thus, 
o-r,vu > 0). With the help of (|2.13p . this equation may be written in the more geometrical 
form [T5] 

sinp f , - a XiW sinv9 f 

- Tbf + Tbf Oi,vu + Til' : 1" o"bf',ri/ <h,™ = 0, (2.17) 

sin (fb 

or 

-Tbf + Tbf Ctr.i/i/ - Tff COS iff 

cos ip h + a^ vv 

-Tff sin iff : h cr b f, rv a v ,Tis = 0, (2.18) 

sin 99b 

where </>f, iff and i^b are the three angles of contact, respectively measured in f, f and b 
(ipi +<pp +(fb = 2-7r). This shows that the classical capillary Young's equation is strongly modified 
and replaced with the preceding one (as it occurred for the thin plate |16yi7j). In the particular 
case of an undeformable solid, owing to a r<uu = 1, a VjTV = and y?b = tt (lim^-j,,,- c °^ h ^ 1 = 0), 
this equation leads to the classical Young's equation —Tbf + Tbf — Tff cos Vf = 0. 



3 Example of Finite-Displacement Solution 

We have shown in the preceding section that the coherence of the theory is mainly based on 
the validity of Green's formula (|2,3p despite the contact line singularity. This point is justified 
in the present section, with an explicit example of solution concerning the simple case of a 
half-space elastic solid, bounded by a plane, and subjected to a constant normal surface tension 
concentrated on a straight line of its surface. This solution satisfies the above line equations 
(|2.13p ~ (|2.18p . its singularity at the contact line is described, its displacement field remains finite, 
the elastic energy is also finite, and it is shown that Green's formula remains valid at the contact 
line. 
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Figure 3: Half-space elastic solid subjected to a normal force concentrated on a straight line of its surface. 
We present a solution with a finite displacement and the formation of an edge at this contact line. 

3.1 A plane strain problem 

Let the body b be an isotropic elastic solid occupying (in the reference state) the half space 
x > 0, in the orthonormal frame (Ox,Oy,Oz'), with a constant and isotropic surface stress of 
eigenvalue a s on its surface x = 0, and subjected to a constant force (per unit length) a\ parallel 
to Ox, concentrated on the line x = y = (Fig. [3j sign convention: cr\ > if the force is directed 
to the outside of b; there is no gravity: g = 0). Clearly, it is equivalent to consider that the body 
is in contact with a fluid f occupying the region x < and y > (in the reference state), and a 
fluid f the other region x < and y < 0, with 75/ = <ti and isotropic surface stresses with the 
same eigenvalue cr s> bf = o" s ,bf = °s- In the present equilibrium state (after deformation), owing 
to the symmetry of the problem with respect to the plane y = 0, and if the surface energies 
7bf and 7bf are the same function of the surface strain tensor e s (temperature and chemical 
potentials being constant), then the preceding equation (|2.17|) is satisfied (a TtUU = 1, a T)TU = 0, 
Tbf = 7bf and iff = (pp, by symmetry). The other equation f|2. 13[) at the contact line gives here 



0"! = 2<7 S COS if, 



(3.1) 



where (p = ipb/2 = ir — ip$ = ir — iff , which determines the angle 92, i.e., the orientation of the 
vector fbf tangent to the bf side of the surface (see Fig. At the surface of the body, instead 
of applying the complex stress condition (|2.8p . we shall impose a simple displacement condition: 



u x = iTTT ' u v = u z < =0 (3.2) 

\y\ + b 

(a 7^ 0, b > 0), at any point (0, y, z') of the surface. The value of a/b 2 is fixed by (|3.1|) : 

75 = d y u x (y = 0+) = 

cr tan if 

= 7= 2 , withp= (3.3) 

In the following (Sec. 13.21 and I3.3p . we solve the problem with b = 1, in the frame of classical 
plane strain elasticity, i.e., with 



u x , u y functions of (x, y) 
u 7 t = 0. 



(3.4) 



By a change of variables (Sec. I3.4p . this will lead to solutions for any a and b satisfying (|3.3p . 
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3.2 The analytic functions F and G 

In the following, z will denote the complex variable x + iy and u the complex displacement 
u x + iu y (function of the complex variable z). We use the general Kolosov's solution of plane 
strain elasticity 

u(z) = -^-(kF(z) + zF 7 {z) + G{zj) 

Zi LI 

where k = — r— — — (3.5) 
X + fj, 

(A, \i > Lame's coefficients, —3 < k < —1), based on the two analytic functions F and 
G. We then follow Muskhelishvili's method (e.g. [2], tome II, annexe XVI) — adapted to the 

i - c 

present singularity problem — to determine F and G. The mapping C — Y z = ui(C) = ^ ^ from 

C — {—1} onto itself is bijective, analytic and w _1 = ui. It transforms B = {C £ C| |C| < 1} into 
A = {z £ C| > 0}, and U - {-1} into D = {z £ C| 3? z = 0} (where U = {C £ C| |C| = 1})- 
The above displacement condition, with b = 1, on the surface of the body thus means 



i.e., with the variable C 



kF(z ) + z F(z ) + G(z Q ) = f(z ) for z £ D, 

where f(z ) = - — — — , y = 5^, (3.6) 
2/0 + 1 



"(Co) 



fc$(Co) + =^$'(Co) + *(Co) = <KCo) forCoGU (3.7) 
"'(Co) 

(extended to Co = -1), where $(C) = F(w(C)), *(0 = ^MO) and 

1 1 



<MCo) = f(z ) 



\yo\ + i *efeg + l 
1 + Co 



l + Co + ie(l-Co)' 

where e = sign(SCo)> for any Co £ U (since |yo| = — ^osignyo = iezo)- The function 6(() 

— n^r = = ■ - — — — is not analytic, but its restriction to U 

w'(C) 1 + C -2 



(3.8) 



#(Co 



i-i (i + Co) 2 l-C 



2 



l-C 2 

is that of the analytic function x(C) = • Since (|3.7|) may be written as 



fc$(Co) + E(Co) = 0(Co) forCoGU, (3.9) 
in which E = x + ^ is analytic, we propose to define $ by 

kH0 = 7^~ [ T^jdCo + C forC£B (3.10) 
2m J 7 Co - C 

(see [H], ibid.), where C is a constant and 7 the circuit t € [0, 2ir] — > e rf . Since is continuous 
in U, $ is analytic in B. Using the decomposition 

<KCo) _ i + Co 1 + C0 

Co-C (1 + Co + *£(!- Co)) (Co -C) (l-ie)(Co + »e)(Co-0 
-1 1 (l + te)(l + C) 1 



ie + C C + ie 2(ie + C) Co-C 
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and according to 



d(o f d( 7T . 

-i 



+ Co + i J 7- Co - i 2 



Co - C J-y- Co - C 
(C G B; 7+ : i G [0, vr] -> e** and 7- : t G [vr, 2vr] e ft ) and 
dCo 



Co-C 



dCo 


/■ d ^ . 




/■ ^ . 


f d( 


Co-C J 


(3+ Co - C V 


Co — c 7 


/3- Co - C J 


a Co - C 



dCo . [ d( f d( 

TTl — I — — TTl 



Co — C J a Co — C J a Co — C 

dz 1-C 
— = -2 log— — , 

«! 2 1 + C 

where /3+ is the path t 6 [0, vr] -> e^*" - ') followed by the path i € [0, 1] ->• 1 + 2t(, (3~ the path 
t G [vr, 2vr] -» e** followed by the path t G [0, 1] ->• 1 + 2*C, « the path i G [0, 1] — >■ 1 + 2(1 - 
ai the path i G [0, 1] — >• 1 + £ — 2t£, and log defined in C — R„, we finally obtain 

km = 2(ilc 2 ) (1+c+c2 + ^" log Trf )+c forCGB ' (3 - n) 

i.e., with the variable z = w(C) 



= Z7TT 2 J 3 + ^ 2 + -^og^) + C 

4(1 + Z Z ) 7T 

§ + z log z 

7r(l + Z 2 ) ' 



(3.12) 



with C = — |. Since | + zlogz = (z — i)F\{z) ((z + z)-Fi(z), respectively) and = j^-F^O^) 
(^t|-F2(^)) respectively), where i*i and i*2 are analytic in a neighbourhood of % (of — i, respec- 
tively), the function F is analytic in C — R_, and then in A. According to lim z ^o k F(z) = |, 
it may be extended as a continuous function in C — IVL, and then in A = A U D. We may write 

kir F(z) = z log z + g(z), i.e. 
(l + z 2 )g(z) = ^-zHogz, (3.13) 

where g may be continuously extended at 0, and after derivation 

kirF'(z) = logz+l+ g'{z), 
2zg(z) + (1 + z 2 )g'{z) = -3z 2 log z - z 2 , (3.14) 

this last equality showing that g' may be continuously extended at 0. We know that F' is 
analytic in C — R (and then continuous in A — {0}) and ()3.14j) shows that z F'(z) and z F'(z) 
may be extended as continuous functions in C — IVL (and then in A). 
Since F(z ) = F(z ) = F(-z ) for z G D, 

7T + Z fog Z - Zq log(-Zo) 



kF(z ) + kF(z Q ) 



7T(1 + Z 2 ) 
7T + Z (-Z£7r) 1 — Js/ol 



vr(l + z 2 ) l-yl 
1 /(*>) (3.15) 



l + |2/o| 
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(with the notations of (|3,6|) and (|3.8|) ). so that (|3.6|) gives 



G(z ) = f(z )-kF(z )-z F'(z ) 

= kF(z ) + z F'(z ) forz GD. (3.16) 

Owing to this expression, we then define G in C — R^ as 

G(z) = kF(z) + zF'(z), (3.17) 

which is analytic in C — R_ (and then in A) and continuous in C — Rl (and then in A). 

3.3 The solution and its singularity 

According to ()3.5[) and ()3.17|) . our solution u is then 

1 



<z) = -^(KHz) + F{z)) + (z + z)F'(z)), (3.18) 

with F given by (|3.12j) . and is a continuous function in A. Its value, for zq € D, 

„( 20 ) = -±KFW + FS)) = -^- T ^ R (3.19) 

(from (|3.15p ) has the form (|3.2p with 6=1. The displacement u is then finite at z = (i.e., at 
the contact line). 

When the variable £ = lo~ 1 (z) tends to —1, \z\ tends to +oo. Using this variable and the 
expression (|3.1ip (with C = —7), we have 



- 2fi u{z) = k(F(z) + F{z)) + {z + z)F'(z) 

- 2S (^ K ) ) + ( i T i + i T |)(l±«^) 

= «(*•«)) - |±f • s^itjP - C 2 ) - d + C 2 ) - "C^T^). 

which tends to when £ tends to —1 (3>(C) tends to owing to (1 — £ 2 ) log = (1 — £ 2 ) log(l — 
C)-(l-C)(l+C)log(l+C); similarly (1- |C| 2 ) log ^ = (1-|C| 2 ) log(l-C)-(l-|C| 2 ) log(l + C) = 
(1 - (CI 2 ) log(l -0-1(1-0(1 + log(l + - i(l - 0(1 + log(l + C) also tends to 0). This 
shows that the displacement u(z) tends to when \z\ tends to +00. 

Kolosov's expressions of the strain and stress tensors components are then obtained from 
(|3.5p and, according to (|3,17D . 

e xx = ^(-(l + k)F'(z)-zF"(z)-G'(z)) 

= 1-^-2(1 + k)F'(z)-(z + z)F"(z)), 

eyy = l-R(-(l + k)F\z) + zF"(z) + G'(z)) 

= ±-®((z + z)F"(z)), 

e xy = ±-^zF"(z) + G'(z)) 

= -L^l + k )F'(z) + (z + z)F"(z)), (3.20) 



11 



a xx = ^(2F f (z)-zF"(z)-G'(z)) 

= ®{(l-k)F'(z)-(z + z)F"(z)), 

o yy = M(2F'(z) + zF"(z) + G'(z)) 

= Vl((3 + k)F'{z) + {z + z)F"(z)), 

a xy = <Z(zF"(z) + G'(z)) 

= %((l + k)F'(z) + (z + z)F"(z)), 
2X 



3?(F'(z)) = (3 + k)$l(F'(z)). (3.21) 



By derivation of (|3.14p . 



knF"(z) = - + g"(z), 

z 

2g(z)+Azg'{z) + {l + z 2 )g"{z) = -6zlogz-5z, (3.22) 

the last equality showing that g" may be continuously extended at 0. The function F" is analytic 
in C — R_ (and then continuous in A — {0}) and (|3.22|) shows that z F"(z) may be extended as 
a continuous function in C — (and then in A) and that 

1 + e~ 2ie 

lim (z + z)F"(z) = — j- , (3.23) 

z— >0, 9 constant K7T 

where 9 = argz. Let tq be > and Vo = {z € C| 3?^ > and \z\ < ro}. The function zF"(z) 
is continuous and then bounded in the compact set Vo, so that (z + z)F"(z) is also bounded in 
Vo — {0}. In addition, the first equality (|3.14|) . where g' is continuous in A and then bounded 
in Vo, and </(0) = (consequence of the second equality (|3,14p ) lead to 

lim^o ^(F'(z)) = +oo 

q 

hm^o, 9 constant $H(F'(z)) = — 

\^(F'(z))\ <ci|logr| + di inVo"-{0} 

%(F'(z)) is bounded in Vo" - {0} (3.24) 

(r = \z\, 9 = argz; c%, d\ constants > 0). 

The expressions (|3.20p - (|3.2ip and the preceding results show that all the components of the 
strain and stress tensors are continuous in A — {0}, bounded in Vo — {0} excepted 

|£xa:|> \<7xx\, \&yy\ and \a z > z >\ < c|logr| + d inV -{0} (3.25) 

(different constants c, d > for each strain or stress tensor component), and 

lim^o £xx = +oo 

1 1 + cos 29 

lim z _j.o, 9 constant &yy 



h m 2— >-0, 9 constant &xy 



2fl klT 

1 (1 + k)9 -sin 20 



2n kn 
lim^^o Vxx, Oyy and a z i z i = +oo 

_ (l + A;)0-sm20 

hm z ^0, 9 constant O xy — T • (3.26) 

/C7T 
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Explicit expressions of F' and F" to be used in (|3.2U|) - (|3.21|) are obtained from (|3.12p : 

1 — irz + z 2 + (1 — z 2 ) log z 



kirF'(z) 



2\2 



l- 7TZ -2z 2 + 3ttz 3 - 3z 4 + (-6z 2 + 2z 4 ) log z 
knF(z) = -^-^ . (3.27) 

3.4 Change of variables 

As noted after (|3.3p . we have solved the problem with the value b = 1, which means (physical) 
dimensionless quantities y = 3?z and z, and, according to (|3.6p . dimensionless quantities F(z) 
and G{z). Physical quantities are obtained by considering the new variable z/b and the new 
functions F(z) = a'F(z/b) and G(z) = a'G(z/b) = kF(z) + z F'(z), where b > is a length and 
a' € R* a force per unit length. With these new functions, the displacement u given by (|3.5|) or 
(I3.18P becomes 

u(z) = a'u(^) (3.28) 
b 

and the components of the strain and stress tensors, Ejj and <Jij given by (|3.2U|) - (|3.2ip . become 



(3.29) 





a' 






~b 




a! 


( z 


aij(z) 


b 





For zq G D, the displacement becomes 



u(z ) = a!u(^-) = ■ 1 (3.30) 



(from (|3.19p ). which has the general form (|3.2|) with a = ^j. 

3.5 Finite elastic energy 

The elastic energy per unit volume is 

A 



2 



(^rrrr + £yy) + M^rrrr + + ^ £ xy) (3.31) 



(using or iij). Since e ra and e xy are continuous and bounded in Vo — {0}, e 2 y and e 2 y are 
integrable in Vo (considered as C R 2 ). Owing to the inequality ()3.25p . e xx and e 2 x are also 
integrable in Vo, which finally implies that (|3.3ip is integrable in Vq. The elastic energy in Vq 



E e \= f {\{e xx + e y y) 2 + ii{e 2 xx + e 2 yy + 2e 2 xy ))dxdy (3.32) 



(per unit length along the normal Oz' to the xy plane) is then well defined and finite. 



3.6 Validity of Green's formula 

The assumption made in Sec. [2] that Green's formula (|2.3j) remains valid, in order to obtain the 
equilibrium equations at the triple contact line, will be now justified using our present finite- 
displacement solution. Since u is continuous in A x R (considered as C R 3 ; we return here to the 
three-dimensional space, according to (|3.4|) ). we consider its variation w = 5u as also continuous 
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in A x R and then bounded in V, where V = Vox]0,/o[ (Iq > 0), so that the components 
Wi 6 L°°(V) (c L 2 (V), since V is bounded). The partial derivatives djUi are d x u x = e xx , 
dyUy = E yy ( written in (|3.20p ) and 

d y u x = ^(-(l-k)F'(z) + zF"(z) + G'(z)) 

= ±-%(2kF'(z) + (z + z)F"(z)), 

d x u y = ^{(l-k)F'(z) + zF"{z) + G'{z)) 

= ±-Z(2F'(z) + (z + z)F"(z)) (3.33) 

(obtained from (|3.5p and (|3. 17j) : djUi = if either i or j refers to the third coordinate z') and 
are all continuous in (A — {0}) x R, and bounded in (Vo — {0}) x R excepted d x u x = e xx 
which satisfies the inequality ()3.25j) in (Vo — {0}) x R. We may then consider their variation 
5(djUi) = djWi as also continuous in (A — {0}) x R, and bounded in (Vo — {0}) x R excepted d x w x 
which will satisfy an inequality similar to (|3.25p in (Vo — {0}) x R, so that all the derivatives 
djWi £ L 2 (V), then the components W{ G H (V). 

Similarly, the components of the stress tensor (Tij E L 2 (V), since they are continuous in 
(A — {0}) x R, and either bounded in (Vo — {0}) x R or satisfying the inequality (|3.25|) in 
(Vo — {0}) x R. Since they are the real or imaginary part of a linear combination of F'{z) and 
(z + z)F"{z) (see (|3.2ip ). their partial derivatives d\Oij [d\ = d x or d y ) will have the form 

8ia i3 = 5? or 9 (A* F"{z) + k 2 (z + z)F"'{z)) (3.34) 

(different constants k\, ki for each 

The expression f|3.22[> . where g" is continuous in A and then bounded in Vo, leads to 

\F"(z)\ <f + d 2 inV^-{0} (3.35) 

(c2, d,2 constants > 0). The derivation of ()3.22p gives 

k7TF"'(z) = -1 + </"(*), 

z z 

6g'(z) + 6zg"(z) + (1 + z 2 )g"'(z) = -61ogz-ll, (3.36) 

the last equality showing that zg"'{z) and zg"'{z) may be continuously extended at 0, and then 
considered as continuous in A, and then bounded in Vo- The expression (|3.36p then leads to 

\(z + z)F"'(z)\ <f + d 3 inVo~-{0} (3.37) 

(03, d-3 constants > 0). The expression ()3.34p and the inequalities (|3,35p and ()3.37p show that 
the partial derivatives d[Oij € L X (V). 

Nevertheless, these derivatives d\Oij L 2 (V) (for I = x or y, and ij = xx, yy, xy or z'z'), so 
that aij H 1 ^). Let us take the example of d x a xx = -3?((1 + k)F"(z) + (z + z)F'"{z))\ 

-kird x a xx = jft( 1+ - + h(z)) 

(using (|3.22|) and (|3.36|) ). where h(z) = (1 + k)g"{z) + (z + z)g"'(z) is continuous and bounded 
in Vo- Thus, d x a xx £ 2 (V) because 



1 + k z + z 2 (k cos 8 — cos 38) 



2 



(»K 3")) 
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is not integrable in Vo- 

Since cr^- g" -ff 1 (V), Green's formula (|2.3p cannot be directly applied on V, as noted in 
Sec. [2J Nevertheless, in the following, we will show that Green's formula remains valid in this 
case. The open set V is bounded by the surfaces S = {z G C| z = iy, < y < ro}x]0, lo[, 
S' = {z G C| z = iy, -r < y < 0}x]0,/ [ and S = ({z G C| > and \z\ = r }x]0,/ [) U 
(Vo x {0}) U (Vo x {Iq}). Since the components of a and w belong to C 1 (V e ), where < e < ro 
and V e = {z G C| > and e < \z\ < ro}x]0,Zo[j Green's formula may be applied on V e 
(with w = on X) 

/ tr(<r* • Dw) dv = - [ div(o-*) -wdv- [ {a*-w)-nda, (3.38) 

JVe We JS £ US' e UC e 

in which S e = {z G C| z = iy, e < y < ro}x]0, lo[, S' £ = {z G C| z = iy, — ro < y < —e}x]0, Iq[, 
C £ = {z G C| $Rz > and \z\ = e}x]0, Zo[ and the unit normal vectors n are directed to the 
inside of V E . Since the components of a and Dw belong to L 2 (V), tr(<7* • Dw) G ^ X (V) which 
implies that Jy £ tr(<T* ■ Dw) dv tends to J v tr(cr* • Du;) dv when e — > 0. Since the components 
of div(cr*) belong to i 1 (V) (the partial derivatives d\Oij G L X (V)) and those of w to L°°(V), 
div(<r*) • w G i : (V) which again implies that J Ve div(c*) • wdv tends to Jy div(<7*) • tt;dw when 
e — > 0. Moreover, according to (|3.25)) and the functions logr and (logr) 2 being integrable in 
[0, 1] (their respective primitives, r log r — r and r(log r) 2 — 2r log r + 2r, tend to when r — > 0), 
the components of a belong to £ 2 (S) and L 2 (S'). Since the components of w also belong to 
L 2 (S U S') (they are continuous and bounded in V), (a* ■ w) ■ n G L X (S U S') which shows that 
Jg eUS / e (<7* • w) ■ nda tends to JsuS'^* ' w ) ' n da when e — > 0. Finally, for z and j fixed, the 
inequality 

| / <7ij Wi rij da\ < (c| log e| + <i) e 7r e Iq 
(from (|3.25p : e constant, | < e in V) shows that 



lim / (a* ■w)-nda = 0. (3.39) 
The limit of p.38p when e — > is then 

[ tr(a* -Dw)dv = - [ div(a*) -wdv - [ (a*-w)-nda, (3.40) 
Jv Jv Js u s' 

i.e. Green's formula on V. Note that some authors [33] proposed that the volume stresses 
produced a resultant force at the contact line. The result ()3.39p expresses that there is no such 
contribution of the volume stresses at the contact line (see also the comment after (|2.13p ) . 



4 Conclusions 

In this paper, which concerns the mechanical surface properties of a deformable body, the general 
surface and contact line equations are first deduced from a variational formulation (see |18j for 
the physical aspects of the theory), by applying Green's formula both in the whole space and on 
the Riemannian surfaces. Despite the singularity at the triple contact line (due to the action of 
the fluid-fluid surface tension on the body), it is assumed that Green's formula remains valid in 
order to obtain the equations at this line. The explicit example of solution given in Sec.[3]justifies 
this assumption. The equations (|2.8|) - (|2.9p at the surfaces are similar to the Cauchy's equations 
for the volume, but involve a special definition of the 'divergence', the surface stress tensor being 
considered as a linear mapping with values in the whole space, then using the tensorial product 
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of the covariant derivatives on the surface and the whole space. The component normal to 
the surface (|2,12p of the divergence equation ()2.8[) is a generalization of the Laplace's equation 
for a fluid-fluid interface. Similar equations were written, e.g., in El HUES]. Note that our 
thermodynamical approach is valid for any deformable body (and not restricted to elasticity). 
There are two equations at the contact line, which represent: (i) the equilibrium of the line fixed 
on the body (|2.13p (equilibrium of the two surface stresses and the fluid-fluid surface tension); 
(ii) the equilibrium relative to the motion of the line with respect to the body ()2,18p . which leads 
to a strong modification of the classical capillary Young's equation (as in the case of the thin 
plate [16U17] ) . Note that (|2.13|) shows that surface stresses are forces which act on a line fixed on 
the body and that the fluid-fluid surface tension is equilibrated by the two surface stresses (and 
not by the volume stresses of the body). This equation (|2.13p suggests a finite displacement 
and the formation of an edge at the contact line, contrary to the infinite-displacement solution 
obtained from classical elasticity [22023] (in which surface properties are not taken into account). 
Recent experiments seem to support this idea [TU]. The existence of such a finite-displacement 
solution is shown with the explicit example of Sec. [3] satisfying the line equations (|2.13p . (|2.18p . 
This elastic solution, based on the approaches of Kolosov and Muskhelishvili — adapted to the 
present singularity problem — and the theory of analytic functions, leads to a description of the 
singularity at the contact line. While the displacement components are continuous functions, 
their first partial derivatives and the strain tensor components are discontinuous, generally 
having different limits when approaching the contact line under different directions (Sec. 13.31 and 
(|3.33p ). This solution also leads to a finite value of the elastic energy (Sec. 13. 5p . whereas this 
energy is infinite in the classical elastic solution [221123] . Owing to the contact line singularity, 
the stress tensor components do not belong to the Sobolev space i? 1 (V). Although Green's 
formula cannot be directly applied in this case, it is shown in Sec. 13.61 that this formula remains 
valid. This result justifies the theory leading to the line equations ([2.13|) - (j2.18|) . It also proves, 
according to (|3.39p . that there is no force contribution of the volume stresses at the contact line 
(contrary to what was proposed in p3]). The importance of Green's formula and its validity for 
a wider class of functions will be presented in a future paper. 



A Eulerian and Lagrangian Surface Quantities 

Let us denote by xq the position of a point of the body in the reference state, x its position 

in the present state, u = x — xq its displacement between the reference state and the present 

state, x' its position in the varied state, w = x' — x = 5x its displacement between the present 

d d 

state and the varied state, <9o« = tt^ an d d% = tt^- Let Sobf> Sbf and Su f respectively be the 

ox ox 1 

bf dividing surfaces in the reference state, the present state and the varied state, and (x^) and 
(x Q ) arbitrary curvilinear coordinates on So,bf and Sbf, respectively, where Greek indices a, /3, 
7,... belong to {1,2} ((xq) and (x a ) must be clearly distinguished from the three-dimensional 
Cartesian coordinates (xq) and (x l ), respectively). The geometrical transformations such that 
Fq : xq — > x, defined in the part of E occupied by the body b, will now be restricted to the bf 
surfaces. We thus have the mappings Fq : xq — > x from So,bf to Sbf, F : x — > x' from Sbf to S' bf , 
jo : xq — > xq from So,bf to E, j : x — > x from Sbf to E, Gq : xq ^ u from So,bf to E and G : x — >■ w 
from Sbf to E, and their respective tangent linear mappings 4>q : cLxq — > dx from T Xo (So,bf) to 
T x (Sbf), 4> ■ dx — > dx' from T x (Sbf) to T,j,/(S bf ), to : dxo — > dxo from Tj; (So,bf) to E, l : dx — > dx 
from T a .(Sbf) to E, ipQ : dxo — > du from T Xo (So,bf) to E and tp : dx — > dw from Tj;(Sbf) to E. 
For arbitrary vectors dxo and dyo £ T xo (So,bf)) the Lagrangian surface strain tensor is defined 

by 

e s (dx , dy ) = \{dx ■ dy - dx ■ dyo) 
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(see footnot^l), which gives 



-{dx* ■ dy - dx* ■ dy ) 

-{dxQ ■ (l + fo)* ■ (to + V>o) • dy - dx* ■ l* q ■ l ■ dy ) (A.l) 



e s = ^(^Q^o + tQ-Vo + ^o-^o) (A.2) 



(e s being here considered as an endomorphism of T xo (So,bf)) 5 i-e., using the coordinates 

1 



e 



a/3 = -(do a Ui d px Q + d 0/3 Ui d 0a x + d Qa Ui dopu 1 ) (A.3) 



(as a covariant tensor). For any vectors dx and dy G T x (Sbf), the Eulerian infinitesimal surface 
strain tensor 5e s is defined by 

5e s {dx,dy) = -5{dx ■ dy) 

= ~{5{dx) ■ dy + dx ■ 5{dy)) 

= -{dx* ■ ip* ■ i ■ dy + dx* ■ i* ■ ip ■ dy), (A.4) 



which gives 



5e s = ^*-o + L*-^) (A.5) 



{5e s as an endomorphism of T^Sbf)), i.e., 



fcs, Q/ 3 = \{d a Wi dpx 1 + dpWi d a x % ) (A.6) 



2 

(as a covariant tensor). It is related to the variation 5e s of e s through 

5e s (dx ,dy ) = ^5(dx ■ dy) = 5e s (dx,dy), 

i.e., 

dx^ ■ 5e s ■ dyo = dx* ■ 5e s ■ dy 

= dxQ ■ (f>Q ■ 5e s ■ <po ■ dyo, 

which gives 

5e s = (f)*Q ■ 5e s • O 

(5e s and 5e s as endomorphisms), i.e., 

$e S ,af3 = doaX^ d Q px v 5s S £ V 

fe s ,ai3 = d a x C dpx^Se s ^ v (A 

(as covariant tensors). 



5e s = O u • 5e s • O 1 (A.7) 



e.g., dx* is the linear form associated to dx, ipg '■ E -¥ T^So.bf) is the adjoint of tpo. 
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Let dxo and dyo be arbitrary vectors of T a;o (So ) bf), dlo and dl the respective lengths of dxQ 
and dx, vq G T xo (So,bf) a unit vector normal to dxo, and f G T x (Sbf) the unit vector normal to 
dx such that ^ • (</>o ■ ^o) > 0. The relation between the areas dao and da of the two parallelograms 
respectively built with (dx$,dyo) and (dx,dy) 



A vq dlo ■ dyo = v dl ■ dy, 

= I rWr • rhn}\ fn 



where ^ = = | det(J • O )| for any isometry I : T^Sbf) -» T XQ (S ,bf), gives 



v dl = A' 1 4>* - vdl. (A.9) 

We then define the Eulerian surface stress tensor a 8 such that the Eulerian surface stress force 
a s - v dl results from the transport by <fto of the Lagrangian surface stress force 7r s • Vq dlo : 

(j)Q -Ws ■ u Q dl = a s ■ vdl, 

which gives, according to (|A.9D . 

7r s = A(j>Q 1 • a s ■ (f>Q U 

a s = A~ l fa ■ tt s • 4>* Q (A. 10) 

(tt s and a s as endomorphisms of T x . (So,bf) and T-^Sbf), respectively), i.e., 

nf = Ad^d^a^ 

af = A- 1 d 0<: x a d 0v xP 4 V (A.11) 

(as contravariant tensors). 

From (fAlU]) and (K7f\i . we have 

tr(7r s -(5e s ) = Ati(cj)Q 1 ■ a s ■ 5e s ■ O ) 
= Ati(a s ■ 5e s ), 

which leads to the Lagrangian and the Eulerian forms of the work of deformation of a surface 
element 

7r s : 5e s dao = a s : 5e s da, (A. 12) 



i.e.. 



K 13 Se s,a/3 da = erf 3 fe s , a/3 da. 

B Surface Equations 

According to ()A.5p . the last term of the equilibrium condition (|2.6|) is first written as 

tr(a s -Ss B ) = tr(^±^. t *-V) 

= tr^-t'.^+trC^^-i*.^). (B.l) 

S being a Riemannian manifold with boundary T, we then apply Green's formula [3] 

divX da = — / X ■ vdl 
s Jr 
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(where X is any vector field of class C 1 on S — which is compact — , v the field of unit vectors on 
T, tangent to S, normal to T and directed to the inside of S, and da and dl are the Riemannian 
measures on S and T, respectively) to the vector field X = a* ■ i* ■ w = al* • w, where 5^ = i ■ cr s , 
if the components of a s and w belong to C 1 (S). At a point x G Sbf, o^* £ T x (Sbf) ® E* (linear 
mapping from E to T x (Sbf)), i.e., 5^* is a section of the vector bundle T(Sbf) ® (Sbf x E*) over 
Sbf. In order to decompose the term divX = div(<7^* • w), we first need to define a covariant 
derivative and a divergence for <t^*. 

In a general way, let us define the covariant derivative V on the vector bundle (T(S b f)*)® 9 ® 
T(S bf )® p ® (S bf x E*)® s ® (S bf x E)® r over S bf (for any p, q,r,s>0), as the tensorial product 
of the covariant derivative V on (T(S b f)*)® 9 <S> T(S b f)® p (for the Levi-Civita connection) and 
the usual derivative d on (S bf x E*)® s ® (S bf x E)® r = S bf x ((E*)® s ® E® r ), i.e., by 

V X {U®V) = {V X U)®V + U ® {d x V), (B.2) 

for any sections X of T(S bf ), U of (T(S bf )*)®« ® T(S bf )® p , and V of (S bf x E*)® s <g> (S bf x E)® r 
(this definition may be justified by using local frames of (T(S b f)*)® 9 ®T(S b f)® p and (S bf x E*)® s 
®(S bf x E)® r ). For any section W of (T(S bf )*) $5<? ® T(S bf )® p <8> (S bf x E*)® s ® (S bf x E)® r , the 
covariant differential of W is then defined as the linear mapping X7W : X — > VxW, so that 
VW is a section of (T(S bf )*)® (<?+1) ® T(S bf ) 0p ® (S bf x E*)® s ® (S bf x E)® r . As an example, 
for a section W of T(S bf ) ® (S bf x E*), VW is a section of T(S bf )* ® T(S bf ) ® (S bf x E*) and, 
by contraction of the covariant index relative to T(S b f)* and the contravariant index relative to 
T(S b f), we thus define divVF, which is a section of S b f x E*. With respect to a local chart of 
Sbf (coordinates (x Q ); Greek indices a, f3, 7,... belong to {1,2}), with the associated frames of 
T(Sbf) and T(S b f)*, and to a basis of the vector space E (coordinates (x*); Latin indices i, j,... 
belong to {1, 2, 3}), with the associated dual basis of E*, we may thus write the components 



{VW)% = dpWf + TfaW? (B.3) 

fi T 

/3 7 l 



(div W)i = dgWf + TLw?, (B.4) 



where are the Christoffel's symbols of the Levi-Civita connection on Sbf. 
With this definition, we may now write 

V(oT ® w) = (Vol*) ® w + SJ* ® iff 

(if) being the usual derivative of w), then 

V(5g* ■ w) = (V5T) • w + oT ■ iff 

(contraction of the covariant index deriving from cK* and the contravariant index deriving from 
w) and 

div(5^* • w) = div(5^*) • w + al* : iff (B.5) 

(contraction of the covariant index deriving from V and the contravariant index deriving from 
Og*). Green's formula applied to X = cT s * ■ w may then be written as 



tr(cr* ■ l* ■ ip) da = — div(cr s *) ■ wda — / (<r s * ■ w) ■ vdl 
s Js Jr 

= —J div(5^*) • w da — J (cr s • v) ■ w dl. (B.6) 
The condition (|2.7|) is then obtained, from (jB.lj) and (|B.6p (since w = on T). 
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Let us now consider a s as a contravariant tensor (convention used in (|2.1U|) ). denote a s = l-o s 
(by contraction of the covariant index of i and the first contravariant index of a s ) and write 

V(i <g> cr s ) = (Vi) ®u s + i® (Vcr s ), 

hence 

Vcf s = (Vi) • a s + i ■ (Vff s ) 

(contraction of the covariant index deriving from t and the first contravariant index deriving 
from <t s ) and 

div<r s = I : a B + l ■ div<r s (B.7) 

(contraction of the covariant index deriving from V and the second contravariant index deriving 
from er s ; divcr s is the usual surface divergence), where 

I = Vt, (B.8) 

i.e., with the components, 

(divo^ = af tp + (diva s ) Q 8 a x\ (B.9) 

with 

(div^) 4 = d^afdaX^ + T^a^dax' 

la/3 = 9 a /3X l — T^p d~ f X % 

(div«r 8 ) a = dpaf + Y%a2 P + T^ar. (B.10) 
Moreover, for any sections X and Y of T(Sbf), we have 

V x {i ®Y) = (VxO ®Y + l® (V X Y), 

hence 

Vx(l-Y) = (V x l)-Y + l-(V x Y) 

(contraction of the covariant index deriving from i and the contravariant index deriving from 
Y), i.e. 

d X Y = (V X i)-Y + V x Y 

= 1:{X®Y) + V X Y (B.ll) 

(Y and V xY being identified to t ■ Y and t ■ (VxY), respectively), which shows that I is 
the second vectorial fundamental form on Sbf (see [5], (20.12.4)), so that (|B.7p and (|B.11|) 
respectively represent the decomposition of div<j s and dxY into the normal component (I : a s 
and I : (X ® Y), respectively) and the tangential component (divcr s and V xY , respectively), 
with respect to S bf . This leads to ([2TTO - ([2TT2l) . 
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